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Cosmic neutrinos provide a unique window into the otherwise-hidden mechanism of particle 
acceleration in astrophysical objects. A flux of high-energy neutrinos was discovered in 2013, 
and the IceCube Collaboration recently reported the likely association of one high-energy 
neutrino with a flare from the relativistic jet of an active galaxy pointed towards the Earth. 
However a combined analysis of many similar active galaxies revealed no excess from the 
broader population, leaving the vast majority of the cosmic neutrino flux unexplained. Here 
we present the likely association of a radio-emitting tidal disruption event, AT2019dsg, with 
a second high-energy neutrino. AT2019dsg was identified as part of our systematic search for 
optical counterparts to high-energy neutrinos with the Zwicky Transient Facility. The proba- 
bility of finding any coincident radio-emitting tidal disruption event by chance is 0.5%, while 
the probability of finding one as bright in bolometric energy flux as AT2019dsg is 0.2%. Our 
electromagnetic observations can be explained through a multi-zone model, with radio anal- 
ysis revealing a central engine, embedded in a UV photosphere, that powers an extended 
synchrotron-emitting outflow. This provides an ideal site for PeV neutrino production. As- 
suming that the association is genuine, our observations suggest that tidal disruption events 


with mildly-relativistic outflows contribute to the cosmic neutrino flux. 


On 2019 October 1, the IceCube Neutrino Observatory! reported the detection of a ~0.2 
PeV neutrino, IC191001A, with an estimated 59% probability of being of astrophysical origin 
based solely on reconstructed energy”. Seven hours later, the direction of the incoming neutrino 
was observed by the Zwicky Transient Facility (ZTF? as part of our neutrino follow-up program. 
The data were processed by our multi-messenger pipeline (see Methods), which performs searches 
for extragalactic transients in spatial and temporal coincidence with high-energy neutrinos’, and 


the radio-emitting tidal disruption event AT2019dsg was identified as a candidate neutrino source. 


Tidal Disruption Events (TDEs) are rare transients that occur when stars pass close to su- 
permassive black holes (SMBHs). Studies have suggested that TDEs are sources of high-energy 


neutrinos and ultra-high energy cosmic rays &7 


, in particular the subset of TDEs with relativistic 
particle jets® ® !% 11. Those TDEs with non-thermal emission are considered the most likely to be 
sources of high-energy neutrinos. AT2019dsg was thus quickly identified as a promising candidate 
neutrino source!?. Given that there are typically <2 radio-emitting TDEs in the entire northern 
sky at any one time, we find that in the 80 sq. deg. of sky observed during the eight neutrino 


follow-up campaigns by ZTF up to March 2020, the probability of finding a radio-detected TDE- 


neutrino association by chance is <0.5%. With the second highest bolometric energy flux of all 
seventeen TDEs detected by ZTF, the probability of finding a TDE at least as bright as AT2019dsg 
by chance is just 0.2%. These calculations are valid for any isotropic distribution, and therefore 
quantify the probability that the AT2019dsg-IC191001A association would arise from atmospheric 
backgrounds. Our program targets four neutrino population hypotheses'’, of which the greatest 
sensitivity is for TDEs (see Methods). Thus, although not directly reflected in the calculation, the 
impact of multiple hypothesis tests on these estimates would be modest. While an atmospheric 
origin for the IC191001A-AT2019dsg association cannot be excluded, the improbability of chance 
temporal and spatial coincidence substantially reinforces the independent energy-based evidence 


of an astrophysical origin for IC191001A, and indicates that any atmospheric origin is unlikely. 


AT2019dsg was discovered'* by ZTF on 2019 April 9, and was classified as a TDE on the 
basis of its optical spectrum!>. This spectrum showed a redshift of z = 0.051, implying a lumi- 
nosity distance Dj ~ 230 Mpc assuming a flat cosmology with Q, = 0.7 and Ho = 70 km s"?! 


Mpc~!. The optical/UV continuum of AT2019dsg is well described by a single blackbody photo- 


sphere with a near-constant temperature!" of 10*°9*9-° K and radius of 10!4°9+°-°3 cm. The peak 


luminosity of 1044°4*°- erg s?! is in the top 10% of the 40 known optical TDEs to date!*, and the 
temperature is in the top 5%. The late-time evolution is consistent with the rapid formation of an 
accretion disk!” !8 (see Figure 1), which would be expected on these relatively short timescales for 
disruptions around higher-mass SMBHs. Indeed the total mass of the host galaxy of AT2019dsg 
is in the top 10% of all optical TDE hosts. Assuming 50% of the host mass is in the bulge, we 


estimate!” a black hole mass of ~ 3 x 10’Mo. 


AT2019dsg was also detected in X-rays, beginning 37 days after discovery. Though the first 
X-ray observation indicated a bright source, with a high X-ray to optical ratio of Ly /Lopt ~0.1, 
this X-ray flux faded extremely rapidly, as shown in Figure 1. This rate of decline is unprece- 
dented, with at least a factor of 50 decrease in X-ray flux over a period of 159 days. Similar to the 
optical/UV emission, the observed X-ray spectrum is consistent with thermal emission, but from 
a blackbody of temperature 1059 K (0.072 + 0.005 keV) and, assuming emission from a circular 
disk, a radius ~ 2 x 10'! cm. As for most X-ray-detected TDEs”®:*!:??, the blackbody radius ap- 
pears much smaller than the Schwarzschild radius (Rs ~ 10" cm) inferred from the galaxy scaling 
relation'?. X-ray emission is generally expected to arise close to the Schwarzschild radius. Small 


emitting areas can arise from an edge-on orientation, because the relativistic velocities at the inner 


disk can Doppler boost a large area of the disk out of the X-ray band. Since our observations probe 
close to the Wien tail of the spectrum, a small temperature decrease due to absorption would also 
yield a significantly underestimated blackbody radius and luminosity??. The exponential decrease 
of the flux could be caused by cooling of the newly-formed TDE accretion disk'* or increasing 


X-ray obscuration. 


Radio observations shown in Figure 2 reveal a third distinct spectral component, namely 
synchrotron emission from non-thermal electrons. We model this emission with a conical geometry 
as expected for outflows (e.g., jets or winds) that are launched from—and collimated by—the inner 
parts of flared accretion disks that emit close to the Eddington limit. Given that electrons are 
typically accelerated with much lower efficiency than protons in astrophysical accelerators”, we 
assume that they carry 10% of the energy carried by relativistic protons (€. = 0.1). We further 
assume that the magnetic fields carry 0.1% of the total energy (eg = 107°), as indicated by radio 
observations of other TDEs™ and supernovae*>. We note that the opening angle for the outflow is 
largely unconstrained. For a half-opening angle, ¢, of 30° we find R = 1.5 x 101° cm in our first 
epoch (41 days after discovery), increasing to R = 7 x 10'$ cm shortly after the neutrino detection 
(177 days after discovery). These radii scale” as R œ [1 — cos(&)] 9/19. The implied expansion 
velocity is roughly constant at v/c = R/ c = 0.12 + 0.01 during the first three epochs, with a 
significant (> 30) acceleration to v/c = 0.21 + 0.02 for the last epoch. These are the velocities 
of the synchrotron-emitting region, and thus provide a lower limit to the velocity at the base of the 
outflow. Indeed even the hotspots of relativistic jets from active galaxies that are frustrated by gas 


in their host galaxy are typically observed?’ to have subrelativistic expansion velocities of ~ 0.1c. 


The inferred outflow energy, E, shows a linear increase from 2.5 x 109 erg to 2 x 10°° erg 
(Figure 2), which would not be expected from models of TDE radio emission which involve a 
single injection of energy”: 9. The constant increase of energy implies a constant injection rate 
at the base of the outflow of approximately 2 x 10*%ergs~'. While some scenarios can yield an 
increase in inferred energy from a single energy injection, none of these are consistent with the 
full set of observed properties. First, a single ejection with a range of velocities could explain the 
observed linear increase of energy with time (the slower ejecta arrive later), but is incompatible 
with the increasing velocity. Second, an increase of the efficiency for conversion of Poynting 
luminosity to relativistic particles is unlikely because the target density that is available to establish 


this conversion is decreasing. And finally, an apparent increase of the inferred energy due to an 


increase of solid angle that emits to our line of sight is only expected for relativistic outflows that 
decelerate. Instead, for AT2019dsg, the observations suggest the presence of a central engine that 
yields continuous energy injection through a coupling of accretion power to the radio emission*”, 
with acceleration in the final radio epoch due to a decrease in the slope of the ambient matter 


density profile. 


* Neutrino emission from AT2019dsg. With this strong evidence for three distinct emission zones 
derived purely from multi-wavelength observations, we consider whether this picture is consistent 
with AT2019dsg being the source of the neutrino IC191001A. In particular, neutrino production 
requires protons to be accelerated to sufficiently high energies, and to collide with a suitably abun- 
dant target. The detection of a single high-energy neutrino implies a mean expectation in the range 
0.05 < N), tot < 4.74 at 90% confidence, where N, tot is the cumulative neutrino expectation for all 
TDEs that ZTF has observed, while for an individual object the expectation will be significantly 
lower*!. AT2019dsg emits f,.1 ~ 0.16 of the population bolometric energy flux, and if we take 


this as a proxy for neutrino emission, we would expect 0.008 < N, < 0.76 for this source. 


Radio observations confirm that particle acceleration is indeed occurring, and that this con- 
tinues without decline through to the detection of the neutrino at ~180 days post-discovery. Given 
that neutrinos typically take a fraction np ~ 0.05 of the parent proton energy, our accelerator must 
be capable of accelerating protons to at least 4 PeV. We evaluate the Hillas criterion”? that the pro- 
ton Larmor radius be less than the system size, to determine whether this is possible. We use our 
estimates for conditions in the synchrotron zone at the time of neutrino detection, with B ~ 0.07 G 
and R ~ 7 x 1016 cm for the near-contemporaneous radio epoch. Taking this as a baseline, we find 
a maximum proton energy of ~160 PeV, far in excess of our requirements. The Hillas criterion 
can also be satisfied within the engine that powers the radio-emitting outflow because the product 


BR is not expected to decrease at smaller radii (e.g. B x R"! for a toroidal configuration). 


The target for neutrino production can be either photons (py interactions) or protons (pp 
interactions). For a photon target, neutrino production occurs above an energy determined by the 
mass of the A resonance, ma. For a thermal spectrum, of temperature T, we then find e, ~ 
Tv (MÄ — m2) /4e,] = 0.3 x (T/10°K) 'PeV. For the UV photosphere of the TDE, we find 
€, ~ 0.8 PeV, while for the compact X-ray source, we find €, ~ 0.05 PeV. Both of these values 


are compatible with the observed neutrino, for which there is a typical uncertainty of one energy 


decade*, so either photon field could serve as a target. For the UV photosphere, we find that 
the mean free path for the parent proton of a PeV neutrino (~ 2 x 1013 cm, see Supplementary 
Information; SI) is much smaller than the photosphere radius, so the UV photosphere is indeed 
optically thick. At smaller radii, the X-rays would overtake the UV photons as dominant scattering 


targets. 


In the multi-zone model, shown in Figure 3, the thermal photons provide a guaranteed target 
for pion production. However hadrons could in principle also serve as a target, leading us to 
consider a single-zone scenario in which the protons are accelerated at the same location as the 
synchrotron-emitting electrons, with the neutrino spectrum following the same intrinsic energy 
power law as the protons and electrons. For pp neutrino production, high target densities of np ~ 
1/(0ppR) ~ 108 cm? would be required for efficient production of neutrinos, where opp is the 
proton-proton cross section and R ~ 10!” cm is the characteristic size of the radio region at the 
time of neutrino production. This high density could be provided by the unbound stellar debris, 
although this component moves with a typical maximum velocity”? of 0.05c , and therefore the 
majority of this debris would have to be swept up with the outflow. Alternatively, the density could 
be provided by pre-existing gas, although since this gas orbits in the sphere of influence of the 


black hole, it would be challenging to satisfy the upper limits on pre-disruption accretion. 


To obtain the expected neutrino flux from this source we have to estimate the energy carried 
by protons (Æp) that are accelerated above the energy threshold needed to produce high-energy 
neutrinos. The outflow energy of 2 x 10°° erg that we derived from the radio observations (Figure 2) 
represent a lower bound to the energy that is available for particle acceleration in a central engine. 
Indeed, the total energy budget for a TDE is set by the mass of the disrupted star, with Erpe ~ 
(1/2) 0.1 Mo c? ~ 10% erg for a solar-mass star. We will assume 1% of this total energy budget 
is carried by relativistic protons, E, ~ 10°! erg. The total energy in muon neutrinos would then 
be FE, tot = (1/8)E, ~ 10° erg for efficient optically-thick pion production, after accounting for 
the pion decay yield and subsequent neutrino flavour oscillations. Convolving this implied energy 
E, tot With the effective area, Aeg, of IceCube's high-energy neutrino alert selection**, we estimate 
the expected number of neutrino alerts. Approximating the sharply-peaked py neutrino spectrum 
as a monoenergetic flux anywhere between 0.2 PeV < €, < 1 PeV, we find N, = (Ev tot /év) (Aeg 
/ 4r D?) ~ 0.03. Thus any optically-thick py scenario would be sufficient to produce the neutrino 


under these assumptions. 


In contrast to a peaked py neutrino spectrum, for pp production the neutrinos would instead 
follow a power law. Many of these neutrinos would then fall below the threshold of IceCube’s alert 
selection. The associated gamma rays would however fall within the sensitive range of gamma-ray 
telescopes, so this scenario could be securely identified through a joint neutrino-gamma ray sig- 
nal. While no gamma-ray emission was measured using the Fermi-LAT telescope for AT2019dsg 
(see Methods), gamma-ray Cherenkov telescopes may be sensitive to the expected gamma-ray sig- 
nal, and the corresponding low-energy (TeV) neutrino emission could confirm a hadronic origin. 
Conversely, the high optical depth of the UV photosphere would absorb any gamma rays accompa- 


10 


nying py neutrino emission. Some contribution from gamma-dark sources is required to explain 


the large astrophysical neutrino flux’. 


Given the different neutrino spectrum expectations, a search for accompanying lower-energy 
neutrinos could be used to probe the conditions at the site of proton interaction. IceCube has 
already searched for correlations between a sample of TDEs and a neutrino dataset dominated 
by lower-energy events, and reported that thermal TDEs account for less than 39% of the diffuse 
astrophysical flux under the assumption of standard candles following a power-law spectrum*. 
The detection of a single TDE-neutrino association with our program would imply that at least 
3% of the astrophysical neutrino alerts arise from the TDE population, fully compatible with these 
pre-existing limits (see SI). 

As TDE discovery rates have increased substantially since the previous IceCube analysis!® 36, 
future searches will be able to study neutrino emission from TDEs with much greater sensitivity. 
A measurement of O(~1-10) TeV neutrinos without accompanying gamma rays would indicate 
that neutrino production is occurring in the X-ray photosphere, rather than in the UV photosphere. 
Indeed, such a detection would confirm the presence of a hidden X-ray source in the first place, 
while our electromagnetic observations cannot. Conversely, a lack of complementary low-energy 
neutrinos and gamma rays implies that only UV photons serve as a target. Neutrinos can uniquely 
serve as probes of the inner region of TDEs, using this novel method of extragalactic neutrino 
tomography. Now that a persistent central engine has been revealed in coincidence with a high- 


energy neutrino, we can begin to shed light on the role of TDEs as astrophysical accelerators. 
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Figure 1: Multi-wavelength lightcurve of AT2019dsg. Error bars represent 10 intervals. The 
upper panel, a, shows the optical photometry from ZTF (in green and red), alongside UV obser- 
vations from Swift-UVOT (in blue, purple and pink). The late-time UV observations show an 
apparent plateau which is not captured by a single power-law decay. The dashed pink line illus- 
trates a canonical t~°/3 power law, while the dotted pink line illustrates an exponentially-decaying 
lightcurve. Neither model describes the UV data well. The lower panel, b, shows the integrated 
X-ray energy flux, from observations with Swift-XRT (in black) and XMM-Newton (in blue), in the 
energy range 0.3-10 keV. Arrows indicated 30 upper limits. The vertical dotted line illustrates the 


arrival of IC191001A. 
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Figure 2: Synchrotron analysis of AT2019dsg. Panel (a) shows radio measurements from 
MeerKAT (1.3 GHz), VLA (2-12 GHz), and AMI (15.5 GHz), at four epochs with times listed 
relative to the first optical detection. The coloured lines show samples from the posterior dis- 
tribution of synchrotron spectra fitted to the measurements at each epoch, and the dashed lines 
trace the best-fit parameters for that epoch. The free parameters are the electron power-law index 
(p = 2.9 + 0.1), the host baseline flux density, plus the magnetic field and radius for each epoch. 
Panel (b) shows the energy at each epoch for a conical outflow geometry with an half-opening 
angle of 30°. The dotted line indicates a linear increase of energy. Panel (c) shows the correspond- 
ing radius for each epoch, with a dotted line illustrating a linear increase. Error bars represent 10 


intervals. 
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Figure 3: Diagram of the three emission zones in AT2019dsg. Panel (a) shows the temporal 
evolution of the three emission regions, while the geometry of these same regions is illustrated in 
panel (b). The size of the region responsible for radio emission (in gray), as well as the blackbody 
radius for the UV-emission (in purple) is derived from data. The Schwarzschild radius is plotted 
in black for a black hole mass of 3 x 10’M.. The white lines in panel (a) represent a continuous 


outflow with velocity c. 
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Methods 


ZTF Neutrino Follow-up 


ZTF routinely images the visible Northern Sky once every three nights to a median depth of 


20.5’, as part of a public survey? 3. 


For our neutrino follow-up program, this wide-field ca- 
dence is supplemented by dedicated Target-of-Opportunity (ToO) observations scheduled through 


the GROWTH ToO Marshal”. 


With ZTF, we have followed up eight neutrinos in the period from survey start on 2018 March 
20 to 2020 March 31, out of a total of 31 neutrino alerts published by IceCube (see SI, Table S1). 
From 2019 June 17, IceCube published neutrino alerts with improved selection criteria to provide 
an elevated alert rate**. In addition to 1 of the 12 alerts under the old selection, ZTF followed up 7 
of the 19 alerts published under the new selection. In general, we aim to follow all well-localised 
neutrinos of likely astrophysical origin reported by IceCube which are visible to ZTF and can be 
observed promptly. Those alerts not observed by ZTF are summarised in the SI (Table S2). Of 
those 23 alerts not followed up by ZTF, the primary reasons were proximity to the Sun (8/23), 
alerts with poor localisation and low astrophysical probability (6/23) and alert retraction (4/23). 
For events which were reported with estimates of astrophysical probability, we chose not to follow 
up those that had both low astrophysical probability (< 50%) and large localisation regions (> 10 
sq. deg.). This astrophysical probability was not reported for high-energy starting events (HESE) 
under the old IceCube alert selection, nor for one recent alert, IC200107A, that was identified 


outside of the standard alert criteria*®. 


Each neutrino localisation region can typically be covered by one or two ZTF observation 
fields. Multiple observations are scheduled for each field, with both g and r filters, and a sepa- 
ration of at least 15 minutes between images. These observations typically last for 300 s, with a 
typical limiting magnitude of 21.0”. ToO observations are typically conducted on the first two 
nights following a neutrino alert, before swapping to serendipitous coverage as part of the public 
survey. Following observations, images are processed by IPAC, and alert packets are generated 


for significant detections from difference images”. 


This alert stream of significant detections is then filtered by our follow-up pipeline built 
within the AMPEL framework"', a platform for realtime analysis of multi-messenger astronomy 


data. Our selection is based on an algorithm for identifying extragalactic transients. We search 
g ying g 
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ZTF data both preceding and following the arrival of the neutrino. In order to identify candidate 


counterparts to the neutrino, we apply the following cuts to ToO and survey data: 


We reject likely subtraction artifacts using machine learning classification and morphology 


cuts®. 


We reject moving objects through matches to known nearby solar system objects*’. We 
further reject moving objects by requiring multiple detections for each candidate (i.e, at the 


same location) separated temporally by at least 15 minutes. 


We remove stellar sources by rejecting detections cross-matched™ to objects with measured 
parallax in GAIA DR2 data”, defined as non-zero parallax with a significance of at least 
30. We further reject likely stars with machine learning classifications**, based on sources 
detected by Pan-STARRS1*’, removing those objects with an estimated stellar probability 
greater than 80%. 


We identify likely AGN by cross-matching to the WISE survey and applying IR color cuts*®. 


We reject detections consistent with low-level AGN variability. 


We require that objects lie within the reported 90% error region to ensure spatial coincidence, 
and that they are detected at least once following the neutrino arrival time to ensure temporal 


coincidence. 


These cuts typically yield ~0.2 candidates per square degree of sky. Promising candidates 


are prioritised for spectroscopic classification, to confirm or rule out a possible association with a 


given neutrino. 


AT2019dsg (R.A.[J2000] = 314.26 deg, Dec[J2000] = +14.20 deg) was spatially-coincident 


with the 90% localisation of the neutrino IC191001A? (R.A. = 314.08*$56 deg, Dec = +12.947 159 


deg), 


at a distance of 1.27 deg to the best-fit position. It was also temporally coincident, being 


detected by ZTF in our ToO observations following the neutrino detection. There were additionally 


three candidate supernovae found in the error region of IC191001A, consistent with background 
expectations. AT2019dsg was the first TDE identified by our pipeline, and the first TDE to be 


reported in coincidence with any high-energy neutrino. 
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Probability of Chance Coincidence 


During the first 18 months of survey operations, ZTF identified 17 TDEs'*, distributed over 28000 
deg of observed sky (the ZTF survey footprint, after removing sources with a Galactic latitude 
|b] < 7). Of these TDEs, each was typically detected for ~6 months!*. We thus estimate that 
the density of ZTF-detected TDEs is approximately 2.0 x107* per sq. deg. of sky in the survey 
footprint at any given time. Our follow-up pipeline requires that any candidate be detected by 
ZTF in ToO observations following a neutrino, in order to establish temporal coincidence. We 
assume that our neutrino pipeline does not have a significantly higher selection efficiency than the 
dedicated ZTF program to identify TDEs'*, and thus that the latter provides a reasonable estimate 


on the background rate of TDEs passing our pipeline. 


Those TDEs with radio detections are considered the most promising candidates for neutrino 
production, as the radio emission serves as a tracer for the particle acceleration required in neutrino 
sources. We can consider the fraction of TDEs which would additionally be detected in radio, 
assuming that all could be observed. Among the ZTF sample of confirmed TDEs, we undertook 
radio follow-up observations with the VLA for 6, of which 2 were detected. Taking this implied 
radio-emitting fraction of 33%, we then find a final density of 5.9 x 10~° radio-emitting TDEs per 
sq. deg. of surveyed sky. 


ZTF has followed-up eight neutrinos up to January 2020, and has covered a combined lo- 
calisation region of 81.05 sq. deg (see SI, Table S1). With this sky area, the expected number 
of coincident radio-detected TDEs across all of our neutrino follow-up campaigns is 4.8 x 107°. 
The Poisson probability of observing at least one radio-emitting TDE during our entire neutrino 


follow-up campaign is thus 4.8 x 107°. 


As radio follow-up observations of ZTF TDEs were biased towards those most likely to be 
detectable, this estimate is an overly conservative one. Because the bolometric energy flux de- 
rived from UV/optical observations (i.e., the blackbody luminosity over the square of the distance) 
serves as a proxy the non-thermal emission, TDEs which were bright under this metric were pref- 
erentially selected for radio observations. To avoid this selection bias, we can instead directly use 
this bolometric energy flux as a proxy for neutrino flux to identify the most promising candidates 
for neutrino detection, namely those TDEs which are both nearby and luminous. Of the 17 TDEs 
observed by ZTF, AT2019dsg ranks second in this metric. The probability of finding a TDE in our 


neutrino follow-up campaign with a bolometric energy flux that is at least as high as AT2019dsg is 
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thus 1.9x107?. 


Like most other studies in neutrino astronomy***?, these chance coincidence probability 
estimates do not account for the so-called “look-elsewhere effect” from multiple possible hypothe- 
ses. In our case, the ZTF program has sensitivity to four theoretically-motivated neutrino pop- 
ulation hypotheses (TDEs, core-collapse supernovae, gamma-ray bursts and AGN flares)'5. The 
impact of the testing multiple hypotheses is thus modest, and a chance coincidence explanation for 
AT2019dsg-IC191001A remains unlikely. 


Of these four, it should also be noted that TDEs are the one to which our program is most 
sensitive. As detailed in the SI, follow-up programs are generally most effective in identifying 
neutrino emission from TDEs, since this flux should be dominated by nearby sources which can be 
detected by telescopes such as ZTF. Moreover, given their low rate, any individual TDE-neutrino 
association will be easier to identify than for more abundant populations such as AGN or super- 
novae. TDEs also evolve sufficiently slowly to enable extensive photometric and spectroscopic 
follow-up, in marked contrast to fast transients such as GRB afterglows, leading to a higher detec- 


tion efficiency. 


Optical/UV Observations 


Prior to the detection of IC191001A, AT2019dsg had already been repeatedly detected by ZTF P48 
telescope as part of the public MSIP survey, most recently on 2019 September 28. These data were 
supplemented by photometric observations from the 2m Liverpool Telescope”? and SEDM?!:°? 
photometry” obtained using the P60 telescope on Mt Palomar. ToO observations of the neutrino 
localisation field began on 2019 October 1, 7.4 hours after the neutrino detection. A second set of 
observations were performed the following night. In all of these images AT2019dsg was clearly 


visible. 


UV observations of AT2019dsg were conducted as part of a systematic survey of UV proper- 
ties of all ZTF-identified TDEs**, using the UltraViolet/Optical Telescope” (UVOT) on board the 
Neil Gehrels Swift Observatory (Swift)*. Data were reduced with uvot source using an aperture 
of 7” to capture the entire galaxy (the host flux density was subtracted based on the best-fit galaxy 
model** and uncertainties on this baseline are propagated into the reported UVOT difference pho- 
tometry). The first UV observation was performed 15 days after the optical peak on 2019 May 


17, and a bright source spatially coincident with the TDE was detected. Subsequent observations 
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continued at a cadence of 2-3 days, up to 2019 September 7. In this period, AT2019dsg contin- 
ued to steadily dim. An additional observation occurred shortly before the neutrino detection on 
2019 September 27. Follow-up observations were then triggered by the identification of a possible 
association with IC191001A'?, beginning on 2019 October 5. 


The optical/UV data are summarised in the SI (Table S3). We note that in the final ZTF 
observations, the source appears to redden in the optical bands. This could be a signature of 
reverberation due emission from dust heated by the TDE””5$; this dust can reach a temperature of 
~ 2000 K. An important caveat is that the contrast between the transient emission and the host is 
very small for these late-time optical detections, so the residuals in the difference image may need 
to be corrected to account for small systematic offsets. That can only be investigated when the 
images for this portion of the public survey are published as part of the next ZTF data release. We 
note that the UV observations are not subject to the same uncertainty because even at late times 


the transient UV flux is about an order of magnitude brighter than the host baseline. 


AT2019dsg was first classified as a TDE by ePESSTO+ on on 2019 May 131°, and the red- 
shift of AT2019dsg was measured to be z = 0.051. Further high-resolution spectroscopic obser- 
vations were conducted using the De Veny Spectrograph on the 4.3m Lowell Discovery Telescope 
(LDT, PI: Gezari), the Kast Double Spectrograph on the 3m Lick Observatory Shane Telescope 
(Lick, PI: Foley), and the Low Resolution Imaging Spectrograph on the 10m Keck Telescope 
(Keck, PI: Graham), with the most recent spectrum on 2019 September 25. These spectra con- 
firm that AT2019dsg belongs to the common spectroscopic class of TDEs with Bowen fluorescence 
emission lines and broad Ha emission lines!®. We note that the Ca triplet is also clearly visible in 
our late-time spectra (rest-frame 8498 A, 8542A and 8662 A), so the SMBH mass could in prin- 
ciple be inferred more precisely using higher-resolution spectroscopy of this feature®!. Following 
the identification of AT2019dsg as a candidate neutrino source, additional high-resolution spectra 
of the source were taken with the 200in Hale Telescope Double Spectrograph at Palomar Obser- 
vatory (P200, PI: Kasliwal & Kulkarni) on 2019 October 3 and again with Lick on 2019 October 
5 and 2019 October 29 (shown in Extended Data Figure 1). There is no evidence of any signifi- 
cant spectral evolution between these spectra and the most recent pre-neutrino spectrum from 2019 


September 25, and the spectral evolution of AT2019dsg is consistent with that of other TDEs'*. 
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X-ray Observations 


AT2019dsg was first observed in X-rays on 2019 May 17 by the X-Ray Telescope (XRT)*, also 
on board Swift**, as part of a program to categorise the X-ray properties of TDEs. AT2019dsg 
was detected at high significance at this epoch, with a measured energy flux of Fy ~ 4 x 10? 


2 s! (0.3-10 keV). Observations continued with a cadence of 2-3 days, and indicated 


erg cm 
a sharply-declining X-ray flux (see SI, Table S4). The source was last detected on 2019 June 
14, and not detected again in any of the following observations continuing until 2019 September 
7. An additional observation was performed with the X-ray Multi-Mirror Mission (XMM-Newton) 
telescope on 2019 May 30, in the range 0.3-10 keV. The XMM-Newton EPIC-pn observations (pro- 
grams 082204 and 08425; P.I. Gezari) were taken in Wide window Thin! filter mode and reduced 
using standard techniques with the XMM-Newton® Science Analysis System (SAS). The source 
extraction region was a circle of radius 35 arcsec at the location of the optical transient in the X-ray 
image, and the background was measured using a 108-arcsec circular region (shown in Extended 
Data Figure 2). The XMM spectrum was binned using the GRPPHA command, such that there were 
at least 20 counts contained in each bin. It was then fit (x?/dof = 59.26/65) with the disk black- 
body (diskbb) model with Galactic™ and intrinsic (Ny ~ 4 x 102° cm~?) absorption described 
using the phabs model in XSPEC v12.9.1°. The flux was consistent with those of Swift-XRT, 
and provided a high signal-to-noise X-ray spectrum well-fitted with a single disk temperature of 
Tais = 10°° K (0.072 +0.005 keV), shown in Extended Data Figure 3. Following the identifi- 
cation of AT2019dsg as a candidate counterpart to IC191001A!?, additional X-ray observations 
were triggered. AT2019dsg was again not detected, with the first Swift-XRT observation occurring 
on 2019 October 5. An additional XMM observation on 2019 October 23 yielded a deep upper 
limit of 9 x 10714 erg cm? s?! (0.3-10 keV) using the same thermal model, computed at the 30 
confidence level using the XMM SAS/HEASARC command eregionanalyse. 


Radio Observations 


Four observations of AT2019dsg were performed with the Karl G. Jansky Very Large Array (VLA) 
under project code 19A-395 (PI: van Velzen), on 2019 May 22, June 19, August 8 and October 5. 
The array was in its moderately-extended B configuration (maximum baseline 11 km) for the first 
two epochs, and in its most extended A-configuration (maximum baseline 36 km) for the final two 
epochs. Our first epoch, on May 22, was a detection experiment, and we observed only in the 8— 


12 GHz band. Having established the presence of radio emission, we observed over a broader range 
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of frequencies in the subsequent three epochs, using the 2-4 GHz, 4-8 GHz, and 8—12 GHz bands. 
We used 3C 48 as a bandpass and flux density calibrator on May 22, and 3C 286 for the other three 
epochs. We used the nearby extragalactic sources ICRF J204945.8+1003 14 (at 4-8 and 8—12 GHz) 
and ICRF J203533.9+185705 (at 2-4 GHz) to determine the complex gain solutions, which were 
interpolated to AT2019dsg. We used the VLA pipeline to perform external gain calibration, and 
after removing residual radio frequency interference, we imaged the data within CASA, using 
Briggs weighting with a robust parameter of 1. We split each baseband into multiple frequency 
bins for imaging (1 GHz bins above 4 GHz, and 0.5 GHz bins below that) to provide better sampling 
of the broadband spectrum, allowing more precise constraints on the turnover frequency, and better 


spectral modelling. 


Radio observations of the field of AT2019dsg were also conducted using the AMI Large Ar- 
ray (AMI-LA)® 7, AMI-LA is a radio interferometer comprised of eight 12.8m-diameter antennas 
producing 28 baselines that range from 18m up to 110m, which operates with a 5 GHz bandwidth 
around a central frequency of 15.5 GHz. We observed AT2019dsg on several epochs (see SI, Ta- 
ble S5) for four hours each. Initial data reduction, editing, and calibration of the phase, and flux 
density, was carried out using reduce_dc, a customized AMI data reduction software package®. 
Phase calibration was conducted using short interleaved observations of ICRF J205135.5+174336, 
while for absolute flux density calibration we used 3C 286. Additional flagging and imaging were 
performed using Common Astronomy Software Application (CASA). All of our observations 
showed a source consistent with the location of AT2019dsg. We used the CASA task IMFIT to 


find the source flux and position. 


Further observations of AT2019dsg were conducted with the South African MeerKAT tele- 
scope, on 2019 June 19, July 29, October 5, and November 30, with each session being 2 h 
long. We used ICRF J193925.0-634245 as a flux-density calibrator, and ICRF J213032.8+050217 
as a phase and amplitude calibrator. The initial calibration was done using the IDIA MeerKAT 
pipeline (https://idia-pipelines.github.io/docs/processMeerKAT), which is 
implemented in CASA. The observed band was 860 MHz wide and centred on 1280 MHz. We 
imaged the whole primary beam (~ 1°) using the CLEAN algorithm (CASA: tclean) in order to 


remove sidelobes from the many (unrelated) sources within the primary beam. The total CLEAN 
flux density in the field was 1 Jy, and the peak brightness in the images was about 48 mJy beam * 
(not related to AT2019dsg). Since residual small calibration errors dominated the image rms back- 


ground in the initial images, we self-calibrated the data in both phase and amplitude, with the 
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mean amplitude gain being fixed at unity to minimise any drifting of the flux-density scale. The 
resolution is slightly different in each epoch, but was ~11” north-south, and ~6” east-west. Image 


rms background levels also varied, ranging between 25 and 32 uJy beam * 


. There was no sign 
of extended emission or confusing sources near AT2019dsg. The flux density was determined by 


fitting an elliptical Gaussian with the same geometry as the restoring beam to the images. 


The measured flux densities from our radio observations are reported in Extended Data Fig- 
ure 7. For all radio observations, the reported uncertainties include both the image background rms 


and a 5% fractional calibration uncertainty, added in guadrature. 


Gamma-Ray Observations 


We analysed data from the Fermi Large Area Telescope (Fermi-LAT)” 


, sensitive to gamma rays 
with energies from 20 MeV to greater than 300 GeV. During its sky-survey operations, the pair- 
conversion telescope Fermi-LAT scans the entire sky about every three hours, and can monitor 
the variable gamma-ray sky over short and long timescales. We studied the region of AT2019dsg 
in three different time intervals, motivated by the multi-wavelength behavior of the source. The 
first interval (G1) includes 130 days of observations that include the peak of the optical emis- 
sion from 2019 April 4 to 2019 August 12. The second one (G2) spans 2019 August 12 to 2019 
November 20 and covers the apparent UV plateau and the peak of the radio emission. The third 
period (G3) integrates the whole period between the start of G1 up to 2020 January 31. We use 
the photon event class from Pass 8 Fermi-LAT data (P8R3 SOURCE), and select a 15° x 15° 
Region of Interest (RoI) centered on the AT2019dsg position derived from optical observations, 
with photon energies from 100 MeV to 800 GeV. We use the corresponding LAT instrument re- 
sponse functions P8R3 SOURCE V2 with the recommended spectral models g/l iem v07.fits and 
iso P8R3 SOURCE V2 v1.txt for the Galactic diffuse and isotropic component respectively. To 
minimise contamination from gamma rays produced in the Earth’s upper atmosphere, we require 
an instrumental zenith angle 6 < 90° for all events, in addition to the standard data quality cuts 


suggested by the Fermi Science Support Center (https://fermi.gsfc.nasa.gov/ssc/ 


data/analysis/scitools/data preparation.html). We perform a likelihood anal- 
ysis, binned spatially with 0.1* resolution and 10 logarithmically-spaced bins per energy decade, 
using the Fermi-LAT ScienceTools package (fermitools v1.0.1) along with the fermipy package 
v0.17.47!. 
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A search was already performed within the 90% error region during both the 1-day and 1- 
month period prior to the arrival of the high-energy neutrino”. No new gamma-ray source was 
identified, and there was no significant (> 50) detection for any source from the fourth Fermi-LAT 
point source catalog (4FGL”’). Here, we specifically test a point-source hypothesis at the position 
of AT2019dsg under the assumption of a power-law spectrum. We find TS = 0 for all intervals, 
where TS is twice the difference in the maximum log £ of an ROI model with and without the 
source, where £ is the likelihood of the data given the model. Upper limits for the energy flux 
(integrated over the whole analysis energy range) have been derived for a power-law spectrum 
(dN/dE x E") with photon power-law index I = 2 and are listed in Extended Data Figure 6, 


along with the respective time intervals. 


In all three time intervals, we detect a new, non-catalogued gamma-ray emitter in the Rol at 
a significance > 50. This source lies just outside the IC191001A 90% error region, as indicated in 
Extended Data Figure 4. The source, which we label Fermi-J2113.8+1120, is likely the gamma-ray 
counterpart of the radio-loud object GB6 J2113+1121, classified as a flat-spectrum radio quasar 
with redshift z = 1.6374. The detection of an unrelated gamma-ray blazar within the neutrino 
uncertainty area is consistent with the background estimation. On average 1.5 4FGL gamma-ray 
blazars are expected in 20 sq.deg. A lightcurve analysis (Extended Data Figure 5) reveals that 
the source was last detected in gamma rays one month prior to the IC191001A detection, and was 
not significantly detected again until another month after the neutrino detection. These results are 
compatible with the findings of the realtime follow-up of the region, which found no evidence for 


emission on timescales of 1-day and 1-month prior the detection of IC191001A ”. 


Such a long apparent lag between gamma-ray emission and neutrino emission is disfavored 


by recent studies on the temporal behavior of hadronic processes in blazars’® 77 


, suggesting that the 
blazar is unlikely to have produced the neutrino. Hence, given the lack of any obvious connection 
between the gamma-ray observations of Fermi-J2113.8+1120 and IC191001A, we do not discuss 


this source any further. 


The HAWC observatory also reported a search for transient gamma-ray emission on short 
timescales in the localisation of IC191001A”, and set a limit for their most significant position 
at 95% confidence of E?AN/dE = 3.51 x 10 19(E/TeV) 93 TeV cm? s71, in the energy range 
300 GeV to 100 TeV, for the period from 2019 September 30 05:46:52 UTC to 2019 October 02 
06:03:29 UTC. We note that this search covered a relatively large region of the sky, and thus had 
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a large associated trial factor. A dedicated search at the position of AT2019dsg would be more 
sensitive, especially one that additionally targeted the longer period over which the central engine 


is active. 


Radio Analysis 


The four radio spectral energy distributions (SEDs) for AT2019dsg can be described by syn- 
chrotron emission from a population of relativistic electrons. We assume that the electrons are 
accelerated into a power-law distribution in energy dN./dy x y~?. We expect that the lowest- 
energy electrons emit their synchrotron radiation below the synchrotron-self absorption frequency 
with negligible free-free absorption, and in this case the shape of the radio SED is determined by 
just 3 free parameters, the electron power-law index p, the magnetic field B and the source radius 
R: 

F, snelt) a MTT (1 e hv RW) (1) 
here 7, and k, are the emission and absorption coefficients, respectively. The normalization of 


Eq. 1 depends on the source geometry and the so-called microphysical parameters (ce, €g) which 
will be treated separately below. Similar to the case of radio-emitting TDE ASASSN-141i/*%, 
we might expect some steady radio emission from the host galaxy. This baseline flux density is 


parameterised as 


F F V Qbaseline 2 
v,baseline — + baseline (Gu) ( ) 
such that the total flux density is given by 

Festal = F baseline + Fose (3) 


The magnetic field and radius are allowed to change for each epoch, while Ftaseline» Ybaseline and p 


are assumed to be constant during our radio observations. 


Using Eq. 3 to describe the synchrotron spectrum, we apply a Markov chain Monte Carlo 
approach to determine a posterior probability distribution of the electron power-law index, as 
well as the peak frequency (eax) and peak flux density (Fear) for each radio epoch (Extended 
Data Figure 8). We allow the measurement variance to be underestimated by some fractional 


amount f (see, e.g., ref. °!), 


The last epoch of VLA observations, which has the best coverage of the optically thin part 
of the radio SED, yields p = 3.0 + 0.15 and we use this as a Gaussian prior when fitting all data 
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simultaneously. We use a flat (uninformative) priors for the other parameters and we cap Qpaseline 
at O and Fyaseline at 0.1 mJy (because the baseline flux density should be optically thin, a < 0, and 
cannot exceed the observed post-TDE radio flux density). For the time-independent parameters we 
find: Fbasetine = 0.09 + 0.01 mJy, baseline = —0.2 + 0.1, p = 2.9 + 0.1, and In f = —3.4. We find 
no significant covariance between the baseline flux density parameters and the peak frequency or 


peak flux density. 


To estimate the radius and energy of the radio-emitting region from the posterior distribution 
Of Freak» Ypeak» and p, we use the scaling relations from Barniol Duran, Nakar, & Piran (2013). 
These relations depend on the electron power-law index and we propagate the uncertainty on p 
into our estimates of Req and Eeq. Additional assumptions for the geometry and the microphysical 
parameters are required. For the geometry of the outflow, our default model is two conical emitting 
regions with half-opening angles  =30 deg, which yield an area covering factor fa = 1—cos¢ = 
0.13 and a volume factor fy = 2/(3tan 9) = 1.15 (here we follow the convention” that f4 = 1 


and fy = 4/3 parameterise a spherical outflow in the Newtonian limit). 


The equipartition energy, Eo is obtained under the assumption that the system contains only 
electrons and magnetic fields (both uniformly-distributed) and that the total energy is minimised for 
Ep = (6/11) Ee. However we expect that protons carry the bulk of the energy and we parameterise 
this energy in protons by €. = Ee/ Ep with E, the total energy in relativistic protons. After this ad- 
justment, E., is increased by (1 + 1/e¢.)!!/(3+?”), Finally, for systems that are not in equipartition, 
the energy is larger by a factor? (11/17) 9/7 + (6/17)7!/", with n = (eg/(1 — €g))/(6/11), 
with eg the fraction of total energy that is carried by the magnetic field. Motivated by observations 
of GRB afterglows$? 5, supernovae” and the relativistic TDE Swift J1644+574, we use €. = 0.1 


and eg = 107°. 


From the eguipartition magnetic field strength inferred from the first epoch of radio obser- 
vation (see Extended Data Figure 8) we estimate that the cooling time of the electrons that emit at 
10 GHz is 500 days. For the last epoch, the field strength has decreased by a factor of 3 and now 
the cooling time is an order of magnitude longer. We can thus conclude that, unless the magnetic 
field energy density is much higher than the eguipartition value (€g/ce > 1), the observed slope 


of the optically thin part of the radio SEDs reflects the intrinsic power-law index of the electrons. 


From our synchrotron analysis we also obtain the number density of relativistic electrons, 


which in turn yields a lower limit to the total particle density in the radio region. This estimate 
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is relevant for the pp scenario of pion production, which requires a target density of at least ~ 


10° cm~? to have sufficient optical depth. For the energy and radius of last radio epoch, which was 


0.1 


obtained a few days after the neutrino detection, we find an electron density of 1094*0! cm”? (see 


Extended Data Figure 8). 


We note that there is no evidence of correlation between the X-ray and radio emission in 
AT2019dsg, in contrast to coupling found for ASASSN-141i*. Such a correlation would only be 
expected if the X-ray luminosity of AT2019dsg served as a tracer of disk power, but the rapid 
observed fading indicates that the observed X-ray emission in AT2019dsg is instead driven either 


by varying obscuration or temperature evolution. 
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Figure ED1: Spectroscopic evolution of AT2019dsg. The public classification spectrum taken 
with the NTT" is plotted in blue, along with later spectra from LDT, Lick, Keck and P200. Each 
is plotted in arbitrary units, offset for display purposes. The Balmer lines are highlighted in cyan, 
the Hell lines in gray, and the Bowen fluorescen@6lines (OII at 3760A, NIII at 4100A and 4640A) 


in black. Telluric lines are marked with +. 
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Figure ED2: X-ray count map from XMM-Newton. The image was taken 50 days after discovery. 
The green circle indicates the source region, while the red circular region was used to measure the 
background. The best-fit position derived from optical observations is spatially-coincident with 


the center of the X-ray source region. 
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Figure ED3: Soft X-ray spectrum of AT2019dsg. Panel (a) shows the photon flux as a function of 
energy in blue, from the XMM Newton spectrum taken on 2019 May 30. The spectrum was fitted 
with an absorbed disk blackbody model, shown in black. Panel (b) shows the ratio of the data to 


the model, with the horizontal orange line indicating unity. Error bars represent lo intervals. 
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Figure ED4: LAT count map of the Region Of Interest (ROT). The map shows the integrated 
search period G3, showing the IC191001A 90% localisation region in orange. The position of 
AT2019dsg is marked by a white star. The neutrino best-fit position is marked with a orange ‘x’. 
Two gamma-ray sources are significantly detected (> 5 o) in the ROI but outside the neutrino 
uncertainty region as marked with white crosses. There is no excess consistent with the position 
of AT2019dsg. 
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Figure EDS: LAT lightcurve for Fermi-J2113.8+1120. The flux is derived in the 0.1-800 GeV 
energy range for the source during the time interval G3, with evenly spaced binning of 28 days. 
Vertical error bars represent 1o intervals, horizontal bars denote bin width. 20 upper limits are 
shown for bins with TS<9. The orange dashed vertical line marks the arrival time of IC-191001A. 
Since this source lies outside the reported 90% error region (see Extended Data Figure 4), and 
given that the LAT lightcurve shows no obvious correlation with the neutrino arrival time, we 


conclude that it is unlikely to be associated with the neutrino. 
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Figure ED6: Gamma-ray energy flux upper-limits for AT2019dsg. The values are derived as- 
suming a point-source with power-law index 1'=2.0 at the position of AT2019dsg, integrated over 


the analysis energy range 0.1-800 GeV. 
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Figure ED7: Peak frequency and peak flux density of the radio observations. The time (At) is 
measured in the observer frame relative to MJD 58582.8, the date of discovery for AT2019dsg. 


equipartition: E-/Eg = 11/6 (no protons) fiducial model: €. = 0.1; eg = 1073 
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Figure ED8: Properties of radio-emitting region. These values are inferred from the synchrotron 
peak flux and peak frequency (see Extended Data Figure 7), where R is the region radius, E is the 
non-thermal energy, B is the magnetic field strength and ne is the density of non-thermal electrons. 


Except for At, all quantities are reported as logio with the uncertainty (68% CL) listed in brackets. 
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1 Discovery and Classification History of AT2019dsg 


AT2019dsg was discovered by ZTF on 2019 April 9, initially named ZTF19aapreis, and reported 
on 2019 April 22 as a likely extragalactic transient'. AT2019dsg was subsequently classified as 
a TDE on 2019 May 13 by ePESSTO+?. Radio emission was tentatively reported on 2019 May 
23 by AMI-LA?, and confirmed on 2019 July 26 by e-MERLIN*. In addition to observations as 
part of a systematic ZTF search for TDEs°, the association with IC191001A prompted additional 


follow-up. 


2 Neutrino Production in AT2019dsg 


The Hillas criterion® for a system of magnetic field strength B and physical radius R can be ex- 


pressed asf: 


E B R 
a = 1600 Z I 
PeV i Gauss R 1016 cm NE 0) 


where Z is the particle charge, 8 ~ 0.2 is the outflow velocity in units of c and E is the maximum 


charged-particle energy. In order for particle acceleration to occur, the timescale reguired for 
particle acceleration must be shorter than the associated particle cooling timescale. Previous work 
has found this condition can be satisfied in TDEs for relevant energies/$, although a detailed 


calculation is beyond the scope of this work. 


These accelerated protons then need sufficient target density. For a photon target, with py 
pion production via the A resonance, we expect that neutrino production will occur above a thresh- 
old: 


EE, ~ 10.16 GeV? (2) 


With this constraint, we can derive the necessary photon energies required for a target to produce 
IC191001A. Taking the reconstructed neutrino energy of ~0.2 PeV directly, we find a threshold 
photon target of Æ, >40 eV. However, these reconstructed neutrino energies typically have upper 


bounds an order of magnitude or more above the central estimate’, so the true neutrino energy 


could be substantially higher. For example, with a true neutrino energy of ~1 PeV, we would 


instead require photons E, >8 eV for pion production. 


During pion production roughly half of the energy will be lost through the neutrino-less 7° 


channel", while for the charged pion channel energy is shared roughly equally among the decay 


products 7* + e* + v 4 Vu + Vu!!. Thus ~3/8 of the pion energy is transferred to neutrinos, 


with a 1:2:0 flavour composition at source. However, across the cosmological baseline travelled, 
neutrino oscillations will lead to a mixed 1:1:1 composition on Earth. The IceCube realtime event 
selection is dominated by muon neutrinos, a channel which will carry no more than ~1/8 of the 


pionic energy. Thus we find: 


E 
Ey = fra (3) 


where f, < 1 is the conversion efficiency of proton energy to pion energy. We can derive the mean 


free path, À, for a proton: 


1 


pylly 


A = 


(4) 


with cross section dp, ~ 5 x 10 cm? and photon number density n}. For a blackbody of 
temperature Tgp ~ 1046 K, the mean free path for the parent proton of a 1 PeV neutrino is 
A ~ 2 x 10 cm. Accounting for the fact that each proton interaction will lead to a typical energy 


reduction of 20%, we then find: 


=022 ) 


fr(a) =1— el (5) 


for path x. Equating x with the radius of the UV photosphere x ~ 10'*°cm, we then find that 
each proton or neutron will typically undergo ~ 10 interactions, which would represent a high 
efficiency f+ ~ 0.9. We caution that this estimate is only approximate, and that detailed numerical 


simulations are required to accurately calculate the pion production efficiency”. 


We then calculate the effective area for a single high-energy neutrino, under the assumption 


of a mono-energetic neutrino spectrum which approximates the expectation for py production. 


The effective area for IceCube varies from 50-200 m? for a 0.2-10 PeV neutrino energy. Below 1 
PeV, this corresponds to a roughly-constant threshold of 6 x 1074 erg cm”? for an expectation of 
one neutrino alert. Given the redshift of AT2019dsg, we find a reguired total energy in neutrinos 
E, =4 x 10°! erg to produce a single neutrino alert. We can thus express the expected number of 


detected neutrinos as: 


N, % 0.03 x (6) 


10°°erg 


This expectation would also be valid for any power-law distribution in the same energy range. 


3 Compatibility with existing constraints 


We can estimate the contribution of TDEs to the diffuse neutrino flux that would be required to 
produce an observation of one association with our ZTF follow-up program. As outlined in Ta- 
ble S1, a total of eight neutrino alerts were observed. For all but one of these, IceCube reported 
an estimate of the signalness, i.e., the probability for each to be astrophysical. We note that this 
quantity is not an absolute value, but is rather derived under specific assumptions about the un- 
derlying neutrino source population. Nonetheless, if we take these estimates at face value, and 
assume that the additional event had the reported signalness mean of 0.5, we would expect a total 
of ~4.3 astrophysical neutrinos in our sample. Taking the implied ZTF population expectation of 
0.05 < Nyro < 4.74, we would then require that a fraction 0.01 < f < 1.00 of the astrophysical 
neutrino flux was produced by ZTF-detected TDEs at 90% confidence. 


We can further consider the contribution of those TDEs that ZTF has not detected, to estimate 
the cumulative contribution of all TDEs to the diffuse neutrino flux. The IceCube collaboration 
has already constrained the contribution of such TDEs to be less than 39% of the total, under the 
assumption of an unbroken E~*° power law and a negative source evolution !*!5. We follow the 
same convention, with any power-law contribution of a transient source population to the diffuse 


neutrino flux is given by: 


1.0 


0.8 


0.6 


CDF 


0.4 


0.2 


0.0 
0.0 0.2 0.4 0.6 0.8 1.0 


Redshift (z) 


Figure S1: Cumulative distribution function (CDF) for TDE neutrino emission as a function 
of redshift. These values are derived under the assumption of an E? power law and a negative 


source evolution!?. Roughly half of all neutrinos will come from sources with z < 0.2. 
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where p(z) is the source rate density, do is the time-integrated particle flux normalisation at 
reference energy Eo. We can use this to calculate the cumulative distribution of neutrino flux as 
a function of redshift'*. This CDF is illustrated in Figure S1 for an E~?" power law, though the 


distribution only depends weakly on the assumed neutrino spectrum. 


It is clear that, for this negative source evolution, the vast majority of TDE neutrinos are 
expected to arrive from sources in the local universe. This statement is independent of both the 
overall level of TDE neutrino production and the absolute TDE rate. ZTF has, thus far, reported the 


detection of TDEs up to a maximum redshift 0.2125. If we simply assume that ZTF can routinely 


detect TDEs up to a redshift of 0.15, fully 40% of the total population flux should come from 
ZTF-detected sources. We would thus require that at least 2.8% of the astrophysical neutrino 
alerts are produced by TDEs, which is fully compatible with the previous IceCube limit. TDE- 
neutrino associations can thus be detected even if the vast majority of the astrophysical neutrino 
flux is produced by other source classes. The large relative contribution of detectable TDEs to the 
population neutrino flux is in marked contrast to supernova-like populations which are dominated 
by distant sources!*, so follow-up searches for TDEs are significantly more sensitive than for these 


other potential sources. 


In any case, IceCube limits are derived under the assumptions of unbroken power laws which 
extend across a broad energy range (100 GeV - 10 PeV), where many additional neutrinos would 
be expected at lower energies. However, for the case of neutrino spectra dominated by high-energy 
components (as expected for py neutrino production), no such low-energy neutrinos would be 


expected, and these existing constraints would then be substantially weakened. 


Supplementary References 


1. Nordin, J. et al. ZTF Transient Discovery Report for 2019-04-22. Transient Name Server 
Discovery Report 2019-615 1 (2019). 


2. Nicholl, M. et al. ePESSTO+ classification of optical transients. The Astronomer’s Telegram 
12752 1 (2019). 


3. Sfaradi, I. et al. A possible radio detection of the TDE candidate AT2019DSG by AMI-LA. 
The Astronomer’s Telegram 12798 1 (2019). 


4. Perez-Torres, M. et al. Unambiguous radio detection of the tidal disruption event AT2019dsg 
with e-MERLIN. The Astronomer’s Telegram 12960 1 (2019). 


5. van Velzen, S. et al. Seventeen Tidal Disruption Events from the First Half of ZTF Survey 
Observations: Entering a New Era of Population Studies. arXiv e-prints arXiv:2001.01409 
(2020). 


6. Hillas, A.M. The Origin of Ultra-High-Energy Cosmic Rays. Annu. Rev. Astron. Astrophys. 
22 425-444 (1984). 


10. 


11. 


12. 


13. 


14. 


15. 


16. 


17. 


18. 


19. 


. Senno, N., Murase, K. & Mészáros, P. High-energy Neutrino Flares from X-Ray Bright and 


Dark Tidal Disruption Events. Astrophys. J. 838 3 (2017). 


. Lunardini, C. & Winter, W. High energy neutrinos from the tidal disruption of stars. Phys. 


Rev. D 95 123001 (2017). 


. IceCube Collaboration et al. Multimessenger observations of a flaring blazar coincident with 


high-energy neutrino IceCube-170922A. Science 361 eaat1378 (2018). 


Himmer, S., Riiger, M., Spanier, F. & Winter, W. Simplified Models for Photohadronic 
Interactions in Cosmic Accelerators. Astrophys. J. 721 630—652 (2010). 


Waxman, E. & Bahcall, J. High energy neutrinos from astrophysical sources: An upper bound. 
Phys. Rev. D 59 023002 (1999). 


Stein, R. Search for High-Energy Neutrinos from Populations of Optical Transients. In 36th 
International Cosmic Ray Conference (ICRC2019), volume 36 of International Cosmic Ray 
Conference, 1016 (2019). 


Sun, H., Zhang, B. & Li, Z. Extragalactic High-energy Transients: Event Rate Densities and 
Luminosity Functions. Astrophys. J. 812 33 (2015). 


Stein, R., Necker, J., Bradascio, F. & Garrappa, S. IceCubeOpenSource/flarestack: Titan 
V2.2.3 (2020). URL https://doi.org/10.5281/zenodo.4005800. 


Pan-Starrs Collaboration et al. Search for transient optical counterparts to high-energy IceCube 
neutrinos with Pan-STARRSI. Astron. Astrophys 626 A117 (2019). 


Blaufuss, E. IceCube-190503A - IceCube observation of a high-energy neutrino candidate 
event. GCN Circular 24378 (2019). 


Stein, R. et al. Optical follow-up of IceCube-190503A with ZTF. The Astronomer's Telegram 
12730 1 (2019). 


Blaufuss, E. IceCube-190629A - IceCube observation of a high-energy neutrino candidate 
event. GCN Circular 24910 (2019). 


Stein, R. et al. Optical follow-up of IceCube-190619A with ZTF. The Astronomer's Telegram 
12879 1 (2019). 


20. 


21. 


22; 


23. 


24. 


25. 


26. 


27. 


28. 


29. 


30. 


31. 


32: 


Stein, R. IceCube-190730A - IceCube observation of a high-energy neutrino candidate event. 
GCN Circular 25225 (2019). 


Stein, R. et al. Optical follow-up of IceCube-190730A with ZTF. The Astronomer’s Telegram 
12974 1 (2019). 


Blaufuss, E. IceCube-190922B - IceCube observation of a high-energy neutrino candidate 
event. GCN Circular 25806 (2019). 


Stein, R., Franckowiak, A., Kowalski, M. & Kasliwal, M. A candidate supernova coincident 
with IceCube-190922B from ZTF. The Astronomer’s Telegram 13125 1 (2019). 


Stein, R., Franckowiak, A., Kowalski, M. & Kasliwal, M. IceCube-190922B: Identification of 
a Candidate Supernova from the Zwicky Transient Facility. GCN Circular 25824 (2019). 


Stein, R. IceCube-191001A - IceCube observation of a high-energy neutrino candidate event. 
GCN Circular 25913 (2019). 


Stein, R., Franckowiak, A., Necker, J., Gezari, S. & Velzen, S.v. Candidate Counterparts to 
IceCube-191001A with ZTF. The Astronomer’s Telegram 13160 1 (2019). 


Stein, R., Franckowiak, A., Necker, J. & Suvi Gezari, S.v. IceCube-191001A: Candidate 
Counterparts with the Zwicky Transient Facility. GCN Circular 25929 (2019). 


Stein, R. IceCube-200107A: IceCube observation of a high-energy neutrino candidate event. 
GCN Circular 26655 (2020). 


Stein, R. & Reusch, S. IceCube-200107A: No candidates from the Zwicky Transient Facility. 
GCN Circular 26667 (2020). 


Stein, R. IceCube-200109A: IceCube observation of a high-energy neutrino candidate event. 
GCN Circular 26696 (2020). 


Reusch, S. & Stein, R. IceCube-200109A: Candidate Counterparts from the Zwicky Transient 
Facility. GCN Circular 26747 (2020). 


Lagunas Gualda, C. IceCube-200117A: IceCube observation of a high-energy neutrino candi- 
date event. GCN Circular 26802 (2020). 


33. 


34. 


35. 


36. 


37. 


38. 


39. 


40. 


41. 


42. 


43. 


44. 


45. 


46. 


Reusch, S. & Stein, R. IceCube-200117A: Candidate Counterpart from the Zwicky Transient 
Facility. GCN Circular 26813 (2020). 


Reusch, S. & Stein, R. IceCube-200117A: One Additional Candidate Counterpart from the 
Zwicky Transient Facility. GCN Circular 26816 (2020). 


Kopper, C. Retraction of IceCube GCN/AMON NOTICE 71165249 130949. GCN Circular 
22669 (2018). 


Blaufuss, E. IceCube-181031A retraction. GCN Circular 23398 (2018). 


Blaufuss, E. Retraction of IceCube GCN/AMON NOTICE 36142391 132143. GCN Circular 
23876 (2019). 


Blaufuss, E. IceCube 41485283 132628.amon retraction. GCN Circular 24674 (2019). 


Blaufuss, E. IceCube-180908A - IceCube observation of a high-energy neutrino candidate 
event. GCN Circular 23214 (2018). 


Taboada, I. IceCube-181014A - IceCube observation of a high-energy neutrino candidate 
event. GCN Circular 23338 (2018). 


Blaufuss, E. IceCube-190124A - IceCube observation of a high-energy neutrino candidate 
event. GCN Circular 23785 (2019). 


Santander, M. IceCube-190704A - IceCube observation of a high-energy neutrino candidate 
event. GCN Circular 24981 (2019). 


Blaufuss, E. IceCube-190712A - IceCube observation of a high-energy neutrino candidate 
event. GCN Circular 25057 (2019). 


Santander, M. IceCube-190819A - IceCube observation of a high-energy neutrino candidate 
event. GCN Circular 25402 (2019). 


Blaufuss, E. IceCube-191119A - IceCube observation of a high-energy neutrino candidate 
event. GCN Circular 26258 (2019). 


Stein, R. IceCube-200227A: IceCube observation of a high-energy neutrino candidate event. 
GCN Circular 27235 (2020). 


10 


47. 


48. 


49. 


50. 


51. 


52. 


53. 


54. 


55. 


56. 


57. 


58. 


Stein, R. IceCube-191215A - IceCube observation of a high-energy neutrino candidate event. 
GCN Circular 26435 (2019). 


Kopper, C. IceCube-190331A - IceCube observation of a high-energy neutrino candidate 
event. GCN Circular 24028 (2019). 


Kopper, C. IceCube-190504A - IceCube observation of a high-energy neutrino candidate 
event. GCN Circular 24392 (2019). 


Taboada, I. IceCube-190921A - IceCube observation of a high-energy neutrino candidate 
event. GCN Circular 23918 (2019). 


Blaufuss, E. IceCube-190629A - IceCube observation of a high-energy neutrino candidate 
event. GCN Circular 24910 (2019). 


Stein, R. IceCube-190922A - IceCube observation of a high-energy neutrino candidate event. 
GCN Circular 25802 (2019). 


Blaufuss, E. IceCube-191122A - IceCube observation of a high-energy neutrino candidate 
event. GCN Circular 26276 (2019). 


Stein, R. IceCube-191204A - IceCube observation of a high-energy neutrino candidate event. 
GCN Circular 26341 (2019). 


Santander, M. IceCube-191231A: IceCube observation of a high-energy neutrino candidate 
event. GCN Circular 26620 (2019). 


Lagunas Gualda, C. IceCube-200120A: IceCube observation of a high-energy neutrino candi- 
date event. GCN Circular 26832 (2020). 


Blaufuss, E. IceCube-200120A: Event likely due to background. GCN Circular 26874 (2020). 


Blaufuss, E. IceCube-181023A - IceCube observation of a high-energy neutrino candidate 
event. GCN Circular 23375 (2018). 


11 


Event R.A. (J2000) Dec (J2000) 90% area ZTF obs = Signalness Ref 
(deg) (deg) (sq. deg.) (sq. deg.) 

IC190503A 120.28 +6.35 1.94 1.37 36% et 
IC190619A 343.26 +10.73 27.16 21.57 55% 18:12 
IC190730A 225.79 +10.47 5.41 4.52 67% 20,21 
IC190922B 5.76 -1.57 4.48 4.09 51% 22, 23, 24 
IC191001A 314.08 +12.94 25.53 20.56 59% 25, 26, 27 
IC200107A 148.18 +35.46 7.62 6.22 = 28, 29 
IC200109A 164.49 +11.87 22,52. 20.06 TI R 50,31 
IC200117A 116.24 +29.14 2.86 2.66 38% 437 


Table S1: Summary of the eight neutrino alerts followed up by ZTF. IC191001A is highlighted 
in bold. The 90% area column indicates the region of sky observed at least twice by ZTF, within the 
reported 90% localisation, and accounting for chip gaps. The signalness estimates the probability 
that each neutrino is of astrophysical origin, rather than arising from atmospheric backgrounds. 


One alert, IC200107A, was reported without a signalness estimate. 


Cause Events 
Alert Retraction IC180423A35, IC181031A3*, IC190205A3”, IC190529A38 
Proximity to Sun IC180908A9, IC181014A*, IC190124A?!, IC190704A 
IC190712A%, IC190819A*, IC191119A*, IC200227A% 
Low Altitude IC191215A" 
Southern Sky IC190331A%, IC190504A9 
Poor Signalness & Localisation IC190221A*9, 1C190629A>!, IC190922A*? 
IC191122A*, 1C191204A*4, IC1912314* 
Bad Weather IC200120A**: 57 
Telescope Maintenance IC181023A 


Table S2: Summary of the 23 neutrino alerts that were not followed up by ZTF since survey 
start on 2018 March 20. Of these, 4/23 were retracted, 11/23 were inaccessible to ZTF for 
various reasons, 6/23 were deemed alerts of poor guality, while just 2/23 were alerts that were 


missed although they passed our criteria. 
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Table S3: Photometry for AT2019dsg, measured by Swift- 
UVOT, ZTF, LT (100) and SEDM. The time (Af) is mea- 
sured in the observer frame relative to MJD 58582.8, the date 
of discovery for AT2019dsg. 


At 
34.02 
38.76 
47.32 
67.26 

0.00 
10.45 
18.10 
23.80 
43.70 
49.43 
52.28 
55.17 
58.98 
64.65 
67.50 
71.30 
72:23 
75.03 
76.11 
77.02 
78.92 
81.71 
93:09 
97.01 
103.60 
104.62 


Band 
i 
i 
i 
i 


r 


Flux [mJy] 
0.19 
0.19 
0.17 
0.13 
0.09 
0.15 
0.19 
0.20 
0.19 
0.16 
0.17 
0.15 
0.15 
0.14 
0.13 
0.11 
0.12 
0.11 
0.11 
0.12 
0.11 
0.11 
0.10 
0.10 
0.09 
0.09 


Flux Error [mJy] 
0.02 
0.02 
0.02 
0.01 
0.01 
0.01 
0.02 
0.02 
0.02 
0.01 
0.02 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 


v [10/4 Hz] 
4.23 
4.23 
4.23 
4.23 
4.96 
4.96 
4.96 
4.96 
4.96 
4.96 
4.96 
4.96 
4.96 
4.96 
4.96 
4.96 
4.96 
4.96 
4.96 
4.96 
4.96 
4.96 
4.96 
4.96 
4.96 
4.96 


Instrument 
LT (100) 
LT (100) 
LT (100) 
LT (100) 

ZTF 
ZTF 
ZTF 
ZTF 
ZTF 
ZTF 
ZTF 
ZTF 
ZTF 
ZTF 
ZTF 
ZTF 
ZTF 
ZTF 
ZTF 
ZTF 
ZTF 
ZTF 
ZTF 
ZTF 
ZTF 
ZTF 
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At Band | Flux [mJy] | Flux Error [mJy] | v [10'* Hz] | Instrument 
106.44 r 0.08 0.01 4.96 ZTF 
109.33 r 0.07 0.01 4.96 ZTF 
115.05 r 0.07 0.01 4.96 ZTF 
129.25 r 0.05 0.01 4.96 ZTF 
163.38 r 0.05 0.01 4.96 ZTF 
166.30 r 0.05 0.01 4.96 ZTF 
167.13 r 0.05 0.00 4.96 ZTF 
32.37 r 0.22 0.02 5.14 SEDM 
34.01 r 0.19 0.02 5.14 LT (100) 
38.76 r 0.17 0.02 5.14 LT (100) 
47.32 r 0.16 0.01 5.14 LT (100) 
23.9 | g 0.20 0.02 6.67 ZTF 
33.29 g 0.19 0.02 6.67 ZTF 
43.76 g 0.18 0.02 6.67 ZTF 
49.46 g 0.16 0.01 6.67 ZTF 
49.48 | g 0.17 0.02 6.67 ZTF 
52.32 g 0.15 0.01 6.67 ZTF 
55.16 g 0.15 0.01 6.67 ZTF 
61.83 g 0.15 0.01 6.67 ZTF 
64.68 g 0.12 0.01 6.67 ZTF 
67.48 g 0.11 0.01 6.67 ZTF 
76.06 g 0.10 0.01 6.67 ZTF 
76.09 g 0.10 0.01 6.67 ZTF 
78.95 g 0.10 0.01 6.67 ZTF 
81.79 g 0.09 0.01 6.67 ZTF 
87.48 g 0.08 0.01 6.67 ZTF 
93.21 g 0.08 0.01 6.67 ZTF 
100.70 g 0.08 0.01 6.67 ZTF 
103.56 g 0.08 0.01 6.67 ZTF 
104.59 g 0.07 0.01 6.67 ZTF 
104.59 g 0.07 0.01 6.67 ZTF 
104.60 g 0.07 0.01 6.67 ZTF 
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At Band | Flux [mJy] | Flux Error [mJy] | v [10'* Hz] | Instrument 
104.64 g 0.07 0.01 6.67 ZTF 
106.42 g 0.07 0.01 6.67 ZTF 
120.79 g 0.05 0.01 6.67 ZTF 
123.42 g 0.06 0.01 6.67 ZTF 
135.90 g 0.04 0.00 6.67 ZTF 
142.54 g 0.04 0.01 6.67 ZTF 
156.79 g 0.03 0.01 6.67 ZTF 
159.57 g 0.03 0.01 6.67 ZTF 
163.42 g 0.03 0.00 6.67 ZTF 
166.20 g 0.03 0.00 6.67 ZTF 
166.22 g 0.03 0.00 6.67 ZTF 
167.16 g 0.03 0.00 6.67 ZTF 
168.12 g 0.03 0.00 6.67 ZTF 
34.01 g 0.19 0.02 6.8 LT (100) 
38.76 g 0.19 0.02 6.8 LT (100) 
47.32 g 0.17 0.02 6.8 LT (100) 
67.26 g 0.12 0.01 6.8 LT (100) 
71.07 g 0.11 0.01 6.8 LT (100) 
74.85 g 0.11 0.01 6.8 LT (100) 
35.90 B 0.34 0.04 7.31 Swift-UVOT 
39.59 B 0.24 0.06 7.31 Swift-UVOT 
42.56 B 0.23 0.06 7.31 Swift-UVOT 
45.53 B 0.22 0.06 7.31 Swift-UVOT 
48.60 B 0.19 0.05 7.31 Swift-UVOT 
51.21 B 0.21 0.05 7.31 Swift-UVOT 
54.59 B 0.17 0.04 7.31 Swift-UVOT 
67.65 B 0.26 0.06 7.31 Swift-UVOT 
35.90 U 0.31 0.02 9.18 Swift-UVOT 
39.59 U 0.27 0.04 9.18 Swift-UVOT 
42.56 U 0.29 0.04 9.18 Swift-UVOT 
45.53 U 0.31 0.03 9.18 Swift-UVOT 
48.60 U 0.27 0.03 9.18 Swift-UVOT 
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At Band | Flux [mJy] | Flux Error [mJy] | v [10'4 Hz] | Instrument 
48.69 U 0.23 0.03 9.18 Swift-UVOT 
51.02 U 0.23 0.03 9.18 Swift-UVOT 
51.21 U 0.24 0.03 9.18 Swift-UVOT 
54.59 U 0.23 0.02 9.18 Swift-UVOT 
57.05 U 0.23 0.02 9.18 Swift-UVOT 
59.89 U 0.18 0.02 9.18 Swift-UVOT 
63.10 U 0.24 0.03 9.18 Swift-UVOT 
67.65 U 0.18 0.03 9.18 Swift-UVOT 
67.99 U 0.19 0.02 9.18 Swift-UVOT 
88.18 U 0.15 0.03 9.18 Swift-UVOT 
88.93 U 0.11 0.02 9.18 Swift-UVOT 
93.09 U 0.11 0.02 9.18 Swift-UVOT 
175.15 U 0.06 0.01 9.18 Swift-UVOT 
179.76 U 0.05 0.02 9.18 Swift-UVOT 
184.37 U 0.04 0.02 9.18 Swift-UVOT 
35.90 | UVWI 0.38 0.02 12.6 Swift-UVOT 
39.59 | UVWI 0.36 0.03 12.6 Swift-UVOT 
42.56 | UVWI 0.34 0.02 12.6 Swift-UVOT 
45.52 | UVWI 0.33 0.02 12.6 Swift-UVOT 
48.60 | UVWI 0.34 0.02 12.6 Swift-UVOT 
48.69 | UVWI 0.29 0.02 12.6 Swift-UVOT 
51.02 | UVWI 0.28 0.02 12.6 Swift-UVOT 
51.21 | UVWI 0.30 0.02 12.6 Swift-UVOT 
54.59 | UVWI 0.30 0.02 12.6 Swift-UVOT 
57.05 | UVWI 0.27 0.01 12.6 Swift-UVOT 
59.89 | UVWI 0.26 0.02 12.6 Swift-UVOT 
63.10 | UVWI 0.28 0.02 12.6 Swift-UVOT 
67.65 | UVWI 0.26 0.02 12.6 Swift-UVOT 
67.99 | UVWI 0.20 0.01 12.6 Swift-UVOT 
88.17 | UVWI 0.17 0.02 12.6 Swift-UVOT 
88.93 | UVWI 0.15 0.02 12.6 Swift-UVOT 
93.09 | UVWI 0.15 0.01 12.6 Swift-UVOT 
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At Band | Flux [mJy] | Flux Error [mJy] | v [10'* Hz] | Instrument 
105.73 | UVWI 0.13 0.01 12.6 Swift-UVOT 
107.43 | UVWI 0.11 0.01 12.6 Swift-UVOT 
111.17 | UVWI 0.10 0.01 12.6 Swift-UVOT 
115.15 | UVWI 0.11 0.01 12.6 Swift-UVOT 
118.88 | UVWI 0.12 0.01 12.6 Swift-UVOT 
123.99 | UVWI 0.09 0.01 12.6 Swift-UVOT 
169.90 | UVWI 0.06 0.01 12.6 Swift-UVOT 
175.15 | UVWI 0.06 0.01 12.6 Swift-UVOT 
179.76 | UVWI 0.06 0.01 12.6 Swift-UVOT 
184.37 | UVWI 0.04 0.01 12.6 Swift-UVOT 
194.09 | UVWI 0.07 0.02 12.6 Swift-UVOT 
199.71 | UVWI 0.07 0.01 12.6 Swift-UVOT 
223.47 | UVWI 0.03 0.01 12.6 Swift-UVOT 
329.98 | UVWI 0.03 0.01 12.6 Swift-UVOT 

35.91 | UVM2 0.37 0.01 14.15 Swift-UVOT 
39.59 | UVM2 0.35 0.02 14.15 Swift-UVOT 
42.57 | UVM2 0.35 0.02 14.15 Swift-UVOT 
45.53 | UVM2 0.35 0.02 14.15 Swift-UVOT 
48.60 | UVM2 0.32 0.01 14.15 Swift-UVOT 
48.69 | UVM2 0.33 0.02 14.15 Swift-UVOT 
51.03 | UVM2 0.26 0.01 14.15 Swift-UVOT 
51.22 | UVM2 0.29 0.02 14.15 Swift-UVOT 
54.60 | UVM2 0.28 0.01 14.15 Swift-UVOT 
57.06 | UVM2 0.28 0.01 14.15 Swift-UVOT 
59.90 | UVM2 0.26 0.01 14.15 Swift-UVOT 
63.10 | UVM2 0.25 0.01 14.15 Swift-UVOT 
67.65 | UVM2 0.26 0.02 14.15 Swift-UVOT 
68.00 | UVM2 0.20 0.01 14.15 Swift-UVOT 
88.18 | UVM2 0.15 0.01 14.15 Swift-UVOT 
88.94 | UVM2 0.18 0.01 14.15 Swift-UVOT 
93.10 | UVM2 0.13 0.01 14.15 Swift-UVOT 
105.72 | UVM2 0.12 0.01 14.15 Swift-UVOT 
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At Band | Flux [mJy] | Flux Error [mJy] | v [10'* Hz] | Instrument 
107.43 | UVM2 0.10 0.01 14.15 Swift-UVOT 
111.16 | UVM2 0.10 0.01 14.15 Swift-UVOT 
115.15 | UVM2 0.10 0.01 14.15 Swift-UVOT 
118.87 | UVM2 0.09 0.01 14.15 Swift-UVOT 
123.99 | UVM2 0.07 0.01 14.15 Swift-UVOT 
169.91 | UVM2 0.05 0.01 14.15 Swift-UVOT 
175.16 | UVM2 0.04 0.00 14.15 Swift-UVOT 
179.77 | UVM2 0.05 0.01 14.15 Swift-UVOT 
184.38 | UVM2 0.04 0.01 14.15 Swift-UVOT 
194.09 | UVM2 0.05 0.01 14.15 Swift-UVOT 
199.72 | UVM2 0.04 0.01 14.15 Swift-UVOT 
223.48 | UVM2 0.04 0.01 14.15 Swift-UVOT 
329.97 | UVM2 0.01 0.00 14.15 Swift-UVOT 

35.90 | UVW2 0.52 0.01 15.69 Swift-UVOT 
39.59 | UVW2 0.47 0.03 15.69 Swift-UVOT 
42.56 | UVW2 0.43 0.02 15.69 Swift-UVOT 
45.53 | UVW2 0.41 0.02 15.69 Swift-UVOT 
48.60 | UVW2 0.42 0.02 15.69 Swift-UVOT 
48.69 | UVW2 0.44 0.02 15.69 Swift-UVOT 
51.02 | UVW2 0.40 0.02 15.69 Swift-UVOT 
51.21 | UVW2 0.39 0.02 15.69 Swift-UVOT 
54.59 | UVW2 0.38 0.01 15.69 Swift-UVOT 
57.05 | UVW2 0.36 0.01 15.69 Swift-UVOT 
59.89 | UVW2 0.36 0.01 15.69 Swift-UVOT 
63.10 | UVW2 0.32 0.02 15.69 Swift-UVOT 
67.65 | UVW2 0.29 0.02 15.69 Swift-UVOT 
68.00 | UVW2 0.31 0.01 15.69 Swift-UVOT 
88.18 | UVW2 0.23 0.01 15.69 Swift-UVOT 
88.93 | UVW2 0.21 0.01 15.69 Swift-UVOT 
93.10 | UVW2 0.20 0.01 15.69 Swift-UVOT 
105.72 | UVW2 0.17 0.01 15.69 Swift-UVOT 
107.43 | UVW2 0.13 0.01 15.69 Swift-UVOT 
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At Band | Flux [mJy] | Flux Error [mJy] | v [10'4 Hz] | Instrument 
111.16 | UVW2 0.12 0.01 15.69 Swift-UVOT 
115.15 | UVW2 0.15 0.01 15.69 Swift-UVOT 
118.87 | UVW2 0.13 0.01 15.69 Swift-UVOT 
123.98 | UVW2 0.10 0.01 15.69 Swift-UVOT 
169.90 | UVW2 0.06 0.01 15.69 Swift-UVOT 
175.15 | UVW2 0.06 0.00 15.69 Swift-UVOT 
179.77 | UVW2 0.06 0.01 15.69 Swift-UVOT 
184.38 | UVW2 0.07 0.01 15.69 Swift-UVOT 
194.09 | UVW2 0.06 0.01 15.69 Swift-UVOT 
199.72 | UVW2 0.04 0.01 15.69 Swift-UVOT 
223.47 | UVW2 0.06 0.01 15.69 Swift-UVOT 
329.96 | UVW2 0.02 0.00 15.69 Swift-UVOT 
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Table S4: X-ray observations of AT2019dsg from Swift- 
XRT and XMM-Newton. The time (Af) is measured in the 
observer frame relative to MJD 58582.8. After At = 65.96, 
the source was not detected. For these observations, we in- 


stead report 30 upper limits. 


At Energy Flux Flux Err Energy Range | Instrument 
[10 ? erg cm~? s7t] | [107}? erg cm~? s71] [keV] 

37.37 4.27 0.42 0.3-10 Swift-XRT 
41.24 1.27 0.67 0.3-10 Swift-XRT 
44.37 1.97 0.46 0.3-10 Swift-XRT 
47.48 3.45 0.61 0.3-10 Swift-XRT 
50.16 1.56 0.04 0.3-10 XMM-Newton 
50.75 2.40 0.34 0.3-10 Swift-XRT 
53.36 1.30 0.26 0.3-10 Swift-XRT 
57.01 0.18 0.10 0.3-10 Swift-XRT 

59.6 0.78 0.23 0.3-10 Swift-XRT 
62.59 0.38 0.15 0.3-10 Swift-XRT 
65.96 0.49 0.23 0.3-10 Swift-XRT 
70.94 0.37 - 0.3-10 Swift-XRT 
92.72 0.46 - 0.3-10 Swift-XRT 
97.49 0.34 - 0.3-10 Swift-XRT 
110.76 0.78 - 0.3-10 Swift-XRT 
112.56 0.96 - 0.3-10 Swift-XRT 
116.48 0.79 - 0.3-10 Swift-XRT 
120.67 0.64 - 0.3-10 Swift-XRT 
124.59 0.66 - 0.3-10 Swift-XRT 
129.96 0.84 - 0.3-10 Swift-XRT 
146.44 2.99 - 0.3-10 Swift-XRT 
149.09 0.98 - 0.3-10 Swift-XRT 
150.64 0.81 - 0.3-10 Swift-XRT 
178.23 0.66 - 0.3-10 Swift-XRT 
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At Energy Flux Flux Err Energy Range | Instrument 
[10-2 erg cm~? s~"] | [107}? erg cm~? s71] [keV] 
183.68 0.30 - 0.3-10 Swift-XRT 
196.16 0.09 - 0.3-10 XMM-Newton 
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Table S5: Radio observations of AT2019dsg from MeerKAT, 
VLA, and AMI-LA, grouped into quasi-simultaneous 


epochs. The time (Af) is measured in the observer frame rel- 
ative to MJD 58582.8, the date of discovery for AT2019dsg. 


At | v [GHz] | Flux density [mJy] | Flux Err [mJy] | Instrument 
42 8.49 0.290 0.032 VLA 
42 9.51 0.408 0.035 VLA 
42 10.49 0.440 0.036 VLA 
42 11.51 0.412 0.037 VLA 
41 15.50 0.464 0.045 AMI 
70 1.28 0.104 0.018 MeerKAT 
70 3.50 0.103 0.028 VLA 
70 4.49 0.319 0.040 VLA 
70 5.51 0.324 0.033 VLA 
70 6.49 0.443 0.037 VLA 
70 7.51 0.558 0.041 VLA 
70 8.49 0.680 0.045 VLA 
70 9.51 0.730 0.047 VLA 
70 10.49 0.756 0.048 VLA 
70 11.51 0.771 0.051 VLA 
71 15.50 0.730 0.076 AMI 
111 1.28 0.111 0.017 MeerKAT 
120 2.24 0.249 0.065 VLA 
120 2.76 0.345 0.053 VLA 
120 3.18 0.255 0.049 VLA 
120 3.69 0.419 0.043 VLA 
120 4.74 0.698 0.046 VLA 
120 5.76 0.829 0.053 VLA 
120 6.69 0.987 0.058 VLA 
120 7.71 1.117 0.063 VLA 
120 8.49 1.194 0.067 VLA 
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At | v [GHz] | Flux density [mJy] | Flux Err [mJy] | Instrument 
120 9.51 1.238 0.069 VLA 
120 | 10.14 1.310 0.073 VLA 
120 | 11.16 1.356 0.075 VLA 
119 | 15.50 0.978 0.075 AMI 
179 1.28 0.152 0.019 MeerKAT 
178 2.24 0.351 0.055 VLA 
178 2:15 0.744 0.054 VLA 
178 3.24 0.920 0.057 VLA 
178 3.76 1.032 0.063 VLA 
178 4.74 1.349 0.073 VLA 
178 5.76 1.379 0.075 VLA 
178 6.69 1.285 0.071 VLA 
178 7.71 1.111 0.062 VLA 
178 8.49 1.074 0.062 VLA 
178 9.51 0.921 0.054 VLA 
178 | 10.14 0.884 0.053 VLA 
178 | 11.16 0.785 0.049 VLA 
179 | 15.50 0.676 0.055 AMI 
235 1.28 0.178 0.032 MeerKAT 
236 | 15.50 0.503 0.047 AMI 
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